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ABSTRACT 

The single-degenerate channel for the progenitors of type la supernovae (SNe 
la) are currently accepted, in which a carbon-oxygen white dwarf (CO WD) ac- 
cretes hydrogen-rich material from its companion, increases its mass to the Chan- 
drasekhar mass limit, and then explodes as a SN la. Incorporating the prescription of 
iHachisu et al.l (|l999al) for the accretion efficiency into Eggleton's stellar evolution code 
and assuming that the prescription is valid for all metallicities, we performed binary 
stellar evolution calculations for more than 25,000 close WD binaries with metallici- 
ties Z = 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.004, 0.001, 0.0003 and 0.0001. For our calcu- 
lations, the companions are assumed to be unevolved or slightly evolved stars (WD -|- 
MS). As a result, the initial parameter spaces for SNe la at various Z are presented 
in orbital period-secondary mass (logPi,Af2) planes. Our study shows that both the 
initial mass of the secondary and the initial orbital period increase with metallicity. 
Thus, the minimum mass of the CO WD for SNe la decreases with metallicity Z. The 
difference of the minimum mass may be as large as 0.24 M© for different Z. 

Adopting the results above, we studied the birth rate of SNe la for various Z via 
a binary population synthesis approach. If a single starburst is assumed, SNe la occur 
systemically earlier and the peak value of the birth rate is larger for a high Z. The 
Galactic birth rate from the WD -I- MS channel is lower than (but comparable to) 
that inferred from observations. Our study indicates that supernovae like SN2002ic 
would not occur in extremely low-metallic i ty env ironments, if the delayed dynamical- 
instability model in iHan fc Podsiadlowskil (|2006( ) is appropriate. 
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1 INTRODUCTION 

Type la supernovae (SNe la) play an important role in as- 
trophysics, especially in cosmology. They appear to be good 
cosmological distance indicators and are successfully ap- 
plied to determine cosmological paramet ers (e.g. Q and A; 
iRiess et~aL I ll998l : iPerlmutter et al.lll999D . There is a linear 
relation between the absolute magnitude of SNe la at max- 
imum light and the magnitude drop of B light curve during 
the first 15 days following ma ximum. This relation is now 
known as the Phillips relation (|Phillipsll 19931 ) and is adopted 
when SNe la are used as distance indicators. In this case, 
the Phillips relation is assumed to be valid at high redshift, 
although it comes from low-redshift samples. Obviously, this 
assumption is precarious since the exact nature of SNe la 
is still unclear (see the reviews by iHillebrandt fc Niemeveij 
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I2OOOI : iLeibundgudbOOd '). If the properties of SNe la evolve 
with redshift, the results for cosmology might be different. 
Thus, it is necessary to study the properties of SNe la at high 
redshift. Since metallicity may represent redshift to some ex- 
tent (e.g. metallicity decreases with redshifts), it is a good 
method to study the properties of SNe la at various redshift 
by finding correlations between their properties and metal- 
Ucity. 

It is generally agreed that SNe la originate from the 
thermonuclear runaway of a carbon-oxygen white dwarf (CO 
WD) in a binary system. The CO WD accretes material from 
its companion, increases mass to its maximum stable mass, 
and then explodes as a thermonuclear runaway. Almost half 
of the WD ma ss is converted into radioactive nickel-56 in 
the explosion (|Branch|[2003) . and the amount of nickel- 
56 de termines the maximum luminosity of SNe la (|Arnettl 
Il982h . Some numerical and synthetical results showed that 
metallicity has an effect on the final amount of nickel-56. 
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and thus the maximum luminosity |Timmes et al.l l2003l : 
iTravagho et~aLl l2005l : iPodsiadlowski et al.ll2006l ). There is 
also some other evidence of the correlation between the 
properties of SNe la and m etallicity from observations, e.g . 



Branch fc Berghl[T9 93: Ham uv et al.ll 19961: IWang et al.ll 19971 : 
CappeUaro et al.l]l997i : .Shanks et al.ll2002i . 

Although the fact that SNe la originate from the 
thermonuclear disruption of mass accreting white dwarfs 
is widely accepted, the precise nature of the progeni- 
tor systems remains unclear. According to the nature 
of the companions of the mass accreting white dwarfs, 
two competing scenarios have b een proposed , i.c the 
double-degenerate ch annel (DD, llben fc Tutukovi 11983 : 
IWhelan fc Iben 1987) and the single degenerate channel 



(SD IWhelan fc IbenI [l973l : iNomoto. Thielemann fc Yokoil 
1 1984 1. In the DD channel, two CO WDs with a total mass 
larger than the Chandrasekhar mass limit may coalesce, 
and then explode as a SN la. Although the channel is 
theoretically less favored, e.g. double WD mergers may 
lead to accretion-induced colla pses rather than to SNe la 
ijffiUebrandt fc Niemeveill2000l l. it is premature to exclude 
the channel at present since there exists e vidence that 
the c h annel may con t ribute a few SNe la plowell et al.l 
l2006l : iBranchI l2006l : lOuimbv. Hoflich fc Wheeled l2007l l. 
The single-degenerate Chandrasekh ar model is the most 
widely accepted one at present (IWhelan fc IbenI 1 19731 : 
iNomoto. Thielemann fc Yokoil 1 1984 ). In this model, the 
maximum stable mass of a CO WD is ~ l.y!9>Mp, (close 
to th e Chandrasekhar mass, INomoto. Thielemann fc Yokoil 
1 19841 ). and the companion is probably a main sequence star 



(WD-^RG) (Yungelson et al. 


I995I: Li fc van den Heuvell 


19971: iHachisu et all 1999a 


bl: iNomoto et al.l 119991: 


Langer et all I2OO0I: iHan fc Podsiadlowskil |2004). The 



SD model is supported by some observations. For example, 
variable circumstellar absorpt ion lines were o bserved 
in the spectra of SN la 2006X JPatat et all 120071). which 



indicates the SD nature of its precursor. IPatat et al] (|2007l ') 
suggested that the progenitor of SN 2006X is a WD -I- RG 
system based on the expansion velocity o f the circum- 
stellar material, while iHachisu et al.l 1I2OO8I) argued a WD 
+ MS nature for this SN la. Recentlv. Ivoss fc NelemansI 
l|2008l ) suggested that SN 2007on is also possibly from a 
WD -I- MS channel. Moreover, several WD -I- MS systems 
are possible progenitors of SN la. For example, supersoft 
X-ray sources (SSSs) were suggested as good cand idates for 
the prog enitors of SNe la jHachisu fc Katdl2003al lbh. bome 
of the SS Ss are WD -|- MS systems and some are WD + RG 
systems (|Di Stefano fc Kone' 2003). A direct way to confirm 
the progenitor model is to search for the companion stars 
of SNe la in their remnants. The discovery of the potential 
companion of Tycho's supernova m ay have verified the 
reliability of the WD -f- MS model l|Ruiz-Lapuente et al.l 
I2OO4I : llhara et al.l l2007l ). In this paper, we only focus on 
the WD -I- MS channel, which is a v ery important channel 
for p roducing SNe la in the Galaxy (|Han fc Podsiadlowskil 
I2OO4I ). 

There are two problems for the WD + MS channel when 
confronted with observations, which are the absence of H 
in the vast majority of SNe la spectra and the difficulty 
of producing SNe la with long delay times (older then 10 
Gyr inferred from SNe la in elliptical galaxies in the lo- 



cal universe, iMannucci et al.l |2005|) or extremely sho rt de- 
lay times (shorter than 0.1 Gyr, Mannucci et al.ll200d ). The 
absence of H in the spectra likely implies that the density 
of circumstellar material (CSM) around SNe la is very low 
or the amount of the hydroge n-rich material stripped from 
the c ompanion is very small l|Mattila et al.l [20051 : iLeonardI 
l2007f ). Neither of the two problems can be overcome easily 
in the WD -I- MS channel. We know from observations that 
some SNe la have delay times longer than 10 Gyr or shorter 
than 0.1 Gyr ((Scannapieco fc Bildsten|[2005l : IMannucci et al.l 



I2OO6I : iTotani et al.ll2008l ). but the WD -I- MS channel can 
only account for SNe la wi th delay times shorter than 2 Gyr 
while longer than 0.1 Gyr (|Han fc Podsiadlowskll2004 ). It is 
likely that SNe la with long or short delay times come from 
other channels than the WD -|- MS channel, e.g. SNe la with 
long delay times come from WD -l-RG channel, while those 
with short delay times come from a WD -I- He star channel 
(Wang, Meng, Chen fc Han 2008, in preparation). 

Many works have con c entrat ed o n the WD-I-MS chan- 
nel. IHachisu et all l|l999al lbl l200i ) and INomoto et all (|l999l . 
I200a 7 have studied the WD-I-MS channel by a simple analyt- 
ical method for treating binary i nteractions. S u ch an alytic 
prescriptions, as pointed out by iLanger et all (|2000l ). can 
not describe some mass-transfe r phases, especially those oc - 
curring on a thermal time-scale. iLi fc van den Heuvell (|l997D 
studied this channel from detailed binary evolution calcula- 
tion, while two WD masses, 1.0 and 1.2 Mq, are considered. 
iLanger et al.l (|2000l ) investigated the channel for metallici- 
ties Z = 0.001 and 0.02, but only fo r mass transfer during 
core h ydrogen burning phase (case A). lHan fc Podsiadlowskil 
l|2004l ) carried out a detailed study of the channel including 
case A and early case B (mass transfer occurs at Hertzprung 
gap (HG)) for Z = 0.02. Considering that not all SNe la are 
found in solar metallicity environments, we pay attention to 
the correlation between the properties of SNe la and metal- 
licities here. 

In this paper, we study the WD -I- MS channel com- 
prehensively and systematically at various Z, showing the 
parameter spaces for the progenitors of SNe la and the dis- 
tributions of the initial parameters for the progenitors of SNe 
la. The results can be effectively used in binary population 
synthesis, or be helpful to search for the potential progenitor 
systems of SNe la. In section [21 we simply describe the nu- 
merical code for binary evolution calculations and the grid of 
binary models we have calculated. The evolutionary results 
are shown in section |3] We describe the binary population 
synthesis (BPS) method in section |4] and present the BPS 
results in section [S] In section [S] we briefly discuss our re- 
sults, and finally we summarize the main results in section 

m 



2 BINARY EVOLUTION CALCULATION 



We use the stellar evolution code of lEggletonl (|1971L 
1 19721 . 119731 ) to calculate the binary evolutions of 
WD-I-MS systems. The code has been updated with 
the latest input physics over the la st three decades 
(iHan. Podsiadlowski fc EggletonI 1 1994 IPoIs et all Il995l . 
1 19981 ). Roche lob e overfiow (RLQ F) is treated within the 
code described bv lHan et ai]||2000l ). We set the ratio of mix- 
ing length to local pressure scale height, a — I /Hp, to 2.0, 
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Figure 1. The total accumulation efficiency, r] = ?jHe^Hi ^ ^ 
function of mass-transfer rate, \M2\ - The hydrogen mass fraction, 
X = 0.7, is adopted for the critical accretion rate, Ma- 



and set the convective overshooting pa rameter, (5ov, to 0.12 
l|Pols et al.lll997l : [Schroder et allll997t ). which roughly cor- 
responds to an overshooting length of 0.25-ffp. Ten metallic- 
ities are adopted here (i.e. Z =0.0001, 0.0003, 0.001, 0.004, 
0.01, 0.02, 0.03, 0.04, 0.05 and 0. 06). The opacity tab les for 
these mctallicties are compiled by Chen & Tout ('2007) from 
[igiesias & Rogers (1996) and Alexander & Ferguson ( 199^1). 
For each Z, the initial hydrogen mass fraction is obtained 
by 



X = 0.76 - 3.0Z, 



(1) 



(Pols et al. 1998). Based on the correlation among X, Y and 
Z used here, Pols et al. (1998) well reproduced the color- 
magnitude diagram (CMD) of some clusters. 

Instead of solving stella r structure equations of a WD, 
we adopt the prescription of Hachisu et al.| (|l999ah on WDs 
accreting hydrogen-rich material from their companions. 
The following is a brief introduction of this prescription. 
In the WD -|- MS channel, the companion fills its Roche 
lobe at MS or during HG, and transfers material onto the 
WD. If the mass-transfer rate, IM2I, exceeds a critical value, 
McT, we assume that the accreted hydrogen steadily burns 
on the surface of WD and that the hydrogen-rich material is 
converted into helium at the rate of Mcv The unprocessed 
matter is assumed to be lost from the system as an optically 
thick wind at a rate of Mwind = I M2I — Mrr (IHachisu et ahl 
1 19961 ). Based on the opacity from llglesias fc Rogers! l| 19961 ). 
the optically thick wind is very sensitive to Fe abundance, 
and it is likely that the wind does not work when Z is lower 
than a certain value (i.e. Z < 0.002, jKobavashi et al. 19S3^. 
Thus, there should be an obvious low-metallicity threshold 
for SNe la in comparison with SN II. However, this m etallic- 
ity threshold has not been found (|Prieto et al. I I2OO8I) . Con- 
sidering the uncertainties in the opacities, we therefore as- 
sume rather arbitrarily that the optically thick wind is still 
valid for low metallicities. 

The critical accretion rate is 



where X is hydrogen mass fraction and Mwd is the mass of 
the accreti ng WD (mass is in M p, and mass-accretion rate is 
in Mm /vr. iHachisu et al.llT999al ). We have not included the 
effect of metallicities on equ ation (2) since the effect is very 
small and can be neglected l|Meng et al.ll200^ . 

The following assumptions are adopted when IM2I is 
smaller than Mcr- (1) When IM2I is higher than ^Mcv, the 
hydrogen- shell burning is steady and no mass is lost from 
the system. (2) When \M2\ is lower than but higher 

than |A/cr, a very weak shell flash is triggered but no mass is 
lost from the system. (3) When |Af2| is lower than |Mcr, the 
hydrogen- shell flash is so strong that no material is accumu- 
lated on the surface of the CO WD. We define the growth 
rate of the mass of the helium layer under the hydrogen- 
burning shell as 

Mhc = ?7H|Af2|, (3) 

where ijh is the mass accumulation efficiency for hydrogen 
burning. According to the assumptions above, the values of 
rju are: 

A/„/|A'/2|, |A'/2|>Af„, 

?7H=<^ 1, Afcr ^ lAfat ^ |Afcr, (4) 

0, lAfal < |Afcr. 

Helium is ignited when a certain amount of helium is 
accumulated. If a He- flash occurs, some of the helium is 
blown off from the surface of the CO WD. Then, the mass 
growth rate of the CO WD, AfwD, is 

Mwo = »7HoAfHc = mcValM'z], (5) 

where 77H0 is the mass accumulation effici ency for helium- 
shell flashes, and its value is taken from iKato fc HachisrJ 
(Hooi): 



T/Hc = 1, 

for Mwd < O.SMq- 

-0.35(logA/Hc + 6.1)^ + 1.02, 



log AfHe < -6.34, 
logAfHe ^ -6.34, 



Mcr = 5.3 X 10" 



X 



(A/wD -0.4), 



(2) 



for O.8A^0 ^ AfwD < 0.9Mq. 

( -0.35(log AfHc + 5.6)2 ^ 
?7Hc = < log A/ho < -6.05, 

[ 1, logA;fHe > -6.05, 

for O.9Af0 < Mwo < I.OMq. 

( -0.35(logA/Hc + 5.6)2 + 1.01, 
VKc = < log Afuo < -5.93, 

[ 1, logAfHe > -5.93, 

for I.OMq < AfwD < I.IMq. 

0. 5410gA/Hc + 4.16, 
log Afuo < -5.95, 

VHc = <^ -0.54(log AfHc + 5.6)2 _^ ^ 

-5.95 < logA/Hc < -5.76, 

1, log Afec ^ -5.76, 

for 1.1A^0 < Mwd < I.SMq. 

-0.175(log AfHc + 5.35)2 ^ 
7?Hc = \ log A/ho < -5.83, 

1, log A'/hc > -5.83, 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 
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for 1.3M0 sC AfwD < 1.35Mq. 

{-0.115(logMHc + 5.7)^ + 1.01, 
log A'/ho < -6.05, (12) 
1, log A'/hc ^ -6.05, 

for MwD > 1.35A^0. 

We show the total accumulation efficiency, ry — rjuciju, 
as a function of mass transfer rate in Fig. [T] which clearly 
shows the prescription above. We incorporated this prescrip- 
tion into Eggleton's stellar evolution code and followed the 
evolutions of both the mass donor and the accreting CO 
WD. The mass lost from the system is assumed to take away 
the specific orbital angular momentum of the accreting WD. 
We calculated more than 25,000 WD+MS binary systems 
with various metallicities, and obtained a large, dense model 
grid. The initial masses of donor stars, M2, range from 0.8 
Mq to 4.3 Mq- the initial masses of CO WDs, Af^^, from 
0.629 Mq to 1.20 Mq; the initial orbital periods of binary 
systems, P\ from the minimum value (at which a zero-age 
main-sequence (ZAMS) star fills its Roche lobe) to ~ 20 
days, where the companion star fills its Roche lobe at the end 
of the HG. In the calculations, we assume that the WD ex- 
plodes as a SN la when its m ass reaches the Chandrasekhar 
mass limit, i.e. 1.378 Mq (|Nomoto. Thielemann fc Yokoil 
1 19841 ). 



3 BINARY EVOLUTION RESULTS 

3.1 Initial parameters for the progenitor of SNe la 

To conveniently compare our results with previous studies 
in the literature, we summarize the final outcomes of all the 
binary evolution calculations in the initial orbital period- 
secondary mass (logP', AI2) plane. For the size limit of this 
paper we only show two examples from the ends of the 
metallicity range we explored (i.e. Z = 0.06 and Z = 0.0001) 
in Figs. [5] and O One can download all the figures from 
\bttp ://www.ynao.ac.cn/^ bps/download/xiangcunmeng.htm. 

From the two figures, we see that CO WDs may reach a 
mass of 1.378 Mq during the optically thick wind phase (the 
filled squares) or after optically thick wind while in stable 
(the filled circles) or unstable (the filled triangles) hydrogen- 
burning phase. All these systems are probably progenitors 
of SNe la. Due to strong hydrogen-shell flash or dynami- 
cally unstable mass transfer, many CO WDs fail to reach 
1.378 Mq. Among them, some CO WDs have mass larger 
than 1.30 Mq at the onset of dynamically unstable mass 
transfer (filled stars). These CO WDs still have a chance to 
increase their masses to 1.378 Mq through delayed dynami- 
cal instability (iHan fc Podsiadlowskill200q) and become SNe 
la. These systems are possible progenitors of the SN 2002ic- 
like supernovae (Han & Podsiadlowski 200(t). The delayed 
dynamical instabilit y can successfully explain m any prop- 
erties of SN 2002ic ijHan fc Podsiadlowskil lioO^ , but more 
evidence is necessary to confirm this scenario. Thus, we have 
not included them in the contours of initial parameter space 
for SNe la (FiglH). 

For the contours, the left boundaries are determined by 
the radii of ZAMS stars, i.e. Roche lobe overflow (RLOF) 
starts at zero age, while the systems beyond the right bound- 
aries undergo dynamically unstable mass transfer at the 



base of the red giant branch The upper bound- 

aries are determined by the delayed dynamical instability 
(mainly for Z > 0.004) and the strong hydrogen-shell flash 
(mainly for Z 0.004). When Z > 0.004, the systems above 
the boundaries have mass ratios too large to stabilize mass 
transfer. If Z ^ 0.004, the mass transfer for the systems 
above the boundaries may be dynamically stable, but the 
mass-transfer rate is so high that most of the transferred 
material is lost from the system by the optically thick wind. 
The mass-transfer rate then sharply decreases to less than 
|A/cr after mass ratio inversion. The lower boundaries are 
constrained by the condition that the mass transfer rate is 
larger than ^Mct and that the secondaries have enough ma- 
terial to transfer onto CO WDs, which can then increase 
their masses to 1.378 Mq. 

We show the contours in the (log P\ M2) plane for four 
CO WD masses with various metallicities in Fig. 14] from 
which we may see the influence of metallicity. The results of 
other initial CO WD masses are similar to Fig. [3] We see 
that, with the increasing Z, the contour moves from lower 
left to upper right in the (logP', M2) plane, indicating that 
the progenitor systems have more massive companions and 
longer orbital periods for a higher Z. This is due to the 
correlation between stellar structure and metallicity. Gen- 
erally, high metallicity leads to larger radii of ZAMS stars, 
then the left hand boundary moves to longer initial orbital 
periods. Meanwhile, stars with high metallicity evolve in 
a way similar to those with low metal licity but less mass 
l|Umeda et al.lll999l : [Chen fc Toutll2007^ . Thus, for binaries 
of CO WDs with particular orbital periods, the companion 
mass increases with metallicity. The situation for the right 
hand boundaries is a bit complicated. They move slightly 
to a longer period with Z when Z > 0.04, while the ten- 
dency is reversed when Z < 0.04. This non-monotonic phe- 
nomenon is derived from the effect of composition (deter- 
mined from equation 1) on opacity and, consequently, on 
stellar structure and evolution. For example, with increas- 
ing Z, the radius of a star at the base of RGB decreases when 
Z < 0.04, but increas es when Z ^ 0.04 l|Chen fc Toutll2007l : 
see iMeng et al.ll2008l for details about the effect of equation 
(1) on stellar evolution). 



3.2 The minimum mass of the CO WD for which 
the WD+MS can produce a SN la 

Previous studies show that the minimum mass of CO WDs 
leadi ng to SNe la, Mwn, may be as low as 0.70Af(:) for Z = 
0.02 (|Langer et al.ll2000l : Han fc Podsiadlowskill2004l '). In our 
study, AfwD strongly depends on Z, as shown in Fig. [5] The 
relation between Af^p and log{Z/ZQ) is almost linear, and 
is fltted by 

AfwD/Af© = 0.6679 - 0.08799 log{Z/ZQ), (13) 



^ Systems consisting of a low-mass secondary and a high-mass 
WD possibly undergo dynamically stable mass transfer at the 
base of the RGB and even on the RGB because of the low mass 
ratio. These systems may also contribute to SNe la (WD + RG). 
Since the study of WD -|- RG is beyond the scope of this paper, 
we leave out these systems from the contours. 
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Figure 2. Final outcomes of binary evolution calculation in the initial orbital period-secondary mass (log P'jMj) plane for CO WD 
-|-MS binaries with Z = 0.06, where is the initial orbital period and is the initial mass of donor star (for different initial WD 
masses as indicated in each panel). Filled squares indicate SN la explosions during an optically thick wind phase (IM2I > Mcr). Filled 
circles denote SN la explosions after the wind phase, where hydrogen-shell burning is stable (Mcr ^ \M2\ ^Mcr). Filled triangles 
denote SN la explosions after the wind phase where hydrogen- shell burni ng is mildly unstable (4 Mcr > 1^21 ^ ^^cr). Filled stars 
denote SN la explosions at delayed dynamical instability phase as shown in lHan &: Podsiadlowskl fcOOSi ) . Open circles indicate systems 
that experience novae explosion, preventing the CO WD from reaching 1.378 Mq, while crosses show the systems that are unstable to 
dynamical mass transfer. 



where the relative error of is less than 2.3% in the 

fitting. We see in Fig. [5] that the minimum mass sharply de- 
creases with metallicity and the difference in M™]5 between 
the cases of Z = 0.06 and Z = 0.0001 may be up to 0.24 
Mq . We explain this as follows: As mentioned in subsection 
13.11 for a high Z, the companions in WD -I- MS systems 
for SNe la are more massive, and so have more material to 
transfer onto CO WDs. The CO WDs therefore do not need 
to be as massive for the production of SNe la. Meanwhile, 
the time-scale for mass transfer is the thermal time-scale, 



which increases with metallicity. This leads to a lower mass- 
transfer rate for a high metallicity. As a consequence, less 
material is lost via optically thick wind. This means that a 
larger fraction of material can be accumulated on the sur- 
faces of CO WDs for systems with a high Z. Therefore CO 
WDs can be initially less massive for SNe la for high Z. 
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Figure 3. Similar to Fig E] but for Z = 0.0001. In this fi gure, we do not plot the cases with minimum WD mass for simplicity of 
appearance of the figure. The minimum WD mass is Af^g = 0.860Mq for Z = 0.0001. Out of the region surrounded by the solid line in 
panel (4), there are some filled symbols representing SNe la explosions. We do not include them in the contour since the mass transfer 
for these systems occurs on the giant branch or at the bottom of the giant branch (see footnote 1). 



4 BINARY POPULATION SYNTHESIS 

Adopting the results in section O we have studied the super- 
nova frequency from the WD+M S channel via the rapid bi- 
nary evolution code developed bv lHurlev eTai] ( ,2000 , 2002 ). 
Here after, we use primordial to represent the binaries before 
the formation of WD-(-MS systems and initial for WD-(-MS 
systems. 



4.1 Common envelope in binary evolution 

During binary evolution, the primordial mass ratio (pri- 
mary to secondary) is crucial for the first mass transfer. 
If it is larger than a critical mass ratio, Qc, the first mass 
transfer is d ynamically unsta ble and a common envelope 
(CE) forms ()Paczvhskil 1 19761 ). The ratio t/c varies with 
the evolutiona ry state of the primordial pr i mary at the on- 
set of RLO F ('HicUming & Webbink"l987'; 'Webbink" 19881; 
Han et al.l [200 2: Podsiadlowski et al.. 2002. : ,Chen & Han 
20081 '). In this study, we adopt qc = 4.0 when the primary 
is on MS or HG. This value is supported by detailed binar y 
evolution studies (|Han et al.ll2000l : IChen fc Hanl2002l . 120031 ). 
If the primordial primary is on FGB or AGB, we use 



gc = [1.67~x + 2(ij^)^]/2.13, (14) 

where is the core mass of primordial primary, and x = 
d\n Ri /d\n Mf is the mass- radius exponent of primordial 
primary and varies with composition. If the mass donors 
(primaries) are n aked helium giants , = 0.748 based on 
equation (|14|) (see iHurlev et al.|[200^ for details) . 

Embedded in the CE is a "new" binary consisting of 
the dense core of the primordial primary and the primordial 
secondary. Due to frictional drag within the envelope, the 
orbit of the "new" binary decays and a large part of the 
orbital energy rele ased in the spiral-in process is injected 
into the envelope (|Livio fc Sokeij[l988h . Here, we assume 
that the CE is ejected if 

QCEA_Bo,b = l-Ebindl, (15) 

where AiJorb is the orbital energy released, -Ebind is the 
binding energy of CE, and acE is CE ejection efficiency 
(i.e. the fraction of the released orbital energy used to 
eject the CE). This criterion i s known as standard a - 
formaUsm (|Nelemans et all l2000l : iNelemans fc ToutI 120051 ) . 
Since the internal energy in the envelope is not incorpo- 
rate d into the binding energy, qce may b e greater than 1 
(see lHan. Podsiadlowski fc Eggletonlll995l for details about 
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Figure 4. The contours of initial parameters in (log P'jMj) for different WD masses and different metallicities, in which SNe la are 
expected. The initial masses of WDs are shown in the lower-right region of each panel. Note that, in all the panels, we use the same 
types of lines, as indicated, to represent each metallicity. 



the internal energy). In this paper, we set acE to be 1.0 or 
3.0. 

iNelemans et al.1 (|200GI 1 and iNelemans fc Tout! (|2005l l 
found that it is difficult for the standard a-formalism to 
reproduce a close pair of white dwarfs without the inclu- 
sion of the thermal energy of the envelope. They then sug- 
gested an alternative algorithm equating the angular mo- 
mentum balance (the -7-algorithm. iNele mans & Tout 2005), 
which may ex plain the f ormat ion of all kinds of close bina- 
ries. However, IWebbinkI (j2007h argued that the a-formalism 
with envelope internal energy included can reproduce ob- 
servations best. Note that the formation of WD -I- MS bi- 
naries is not much affected by the use of the 7-algorithm 
in comparison to a-form alism, especially when acE = 3.0 
l|Nelemans fc Toutll2005l ). Therefore, we continue to use the 
a-formahsm throughout this paper and the main results in 
this paper will hold regardless. 



4.2 Evolutionary channels to WD + MS systems 

According to the evolutionary phase of the primordial pri- 
mary at the onset of the first RLOF, WD -f MS systems 
can be produced from three channels. 

Case 1 (He star channel): the primordial primary is in 
HG or on RGB at t he onset of the first RLOF (i. e. case B 
evolution defined bv iKippenhahn fc Weigertlll967^ . In this 
CE is formed because of a large mass ratio or a 
convective envelope of the mass donor. After the CE ejection 
(if it occurs), the mass donor becomes a helium star and 
continues to evolve. The helium star likely fills its Roche 
lobe again after the central helium is exhausted. Since the 
mass donor is much less massive than before, this RLOF is 
dyn amically stable , resul t ing in a close CO WD-I-MS system 
(see iNomoto et al.lll999l . |2003| for details). The primordial 
orbital periods may range from 10 to 1000 days for this 
channel. 
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Figure 5. Relation between the minimum mass of CO WD for 
the production of SNe la and metaUicity, with which the CO WD 
can increase its mass to 1.378 Mq. 



Case 2 (EAGB channel): the primordial primary is in 
early asymptotic giant branch stage (EAGB) (i.e. helium 
is exhausted in the core, while thermal pulses have not yet 
started). A CE is formed because of dynamically unstable 
mass transfer. After the CE is ejected, the orbit decays and 
the primordial primary becomes a helium red giant (HeRG) . 
The HeRG may fill its Roche lobe and start the second 
RLOF. Similar to the He star channel, this RLOF is stable 
and produces WD + MS systems after RLOF. The range of 
the primordial periods is ~ 100 - 1000 days for this channel. 
During the EAGB, the hydrogen-burning shell is extinct in 
the primary, which ha s the same core m ass as that at the 
bottom of the EAGB (|Hurlev et al.ll2000l ). 

Case 3 (TPAGB channel): the primordial primary fills 
its Roche lobe at the thermal pulsing AGB (TPAGB) stage. 
Similar to the above two channels, a CE is formed during 
the RLOF. A CO WD + MS binary is produced after CE 
ejection. The primordial periods of the systems experiencing 
the channel are larger than 1000 days. 

Generally, the CO WD + MS binaries obtained from the 
different channels above have some different properties and 
these properties eventually determine whether the systems 
may explode as SNe la or not. We will discuss this in section 
[5l The WD + MS systems continue to evolve and the MS 
companions may fill their Roche lobes, and transfer their 
material to the CO WDs, which are likely to explode as SNe 
la. Here, we assume that, if the initial orbital period, Porbi 
and the initial secondary mass, M2, of a WD + MS system is 
located in the appropriate regions in the (logP',M2) plane 
(see Figs. [2]to|3 for SNe la at the onset of RLOF, a SN la 
is produced. 



4.3 Basic parameters for Monte Carlo simulations 

To investigate the birth rate of SNe la, we followed the evo- 
lution of lO'^ bi naries for various Z via H urley's rapid binary 
evolution code l|Hurlev et aLll200d . l200j ). The resuhs of grid 
calculations in section|3]are incorporated into the code. The 



code is only valid for Z ^ 0.03, and the production of SNe la 
of seven metallicities, i.e. Z = 0.03, 0.02, 0.01, 0.004, 0.001, 
0.0003 and 0.0001, are investigated here. The primordial bi- 
nary samples are generated in a Monte Carlo way and a 
circular orbit is assumed for all binaries. The basic parame- 
ters for the simulations are as follows. 

(i) The initial mass function (IFM) of iMiUer fc Scald 

^1973) is adopted. The primordial prima r y is generated ac- 
cording to the formula of .Eggleton et al.l l|l989l) 



Mf = 



0.19X 



(l_X)0-75 4_ 0.032(1 -X)° 



(16) 



where X is a random number uniformly distributed in the 
range [0,1] and Mf is the mass of the primordial primary, 
which ranges from 0.1 Mq to 100 Mq. 

(ii) The mass ratio of the primordial binaries, q' , is a 
very important parameter for binary evolution . For simplic- 
ity, we take a uniform mass-ratio distribution l|Mazeh et al.l 
1 19921 : iGoldberg fc Mazeh 1994) : 



n(g') 



< g' ^ 1, 



(17) 



where q' = M^/Mf. 

(iii) We assume that all stars are members of binary 
systems and that the distribution of separations is constant 
in log a for wide binaries, where a is separation, and falls off 
smoothly at small separation: 



a ■ n(a) 



where «= 



(a/aoY 



a < ao; 
ao < a < ai, 



(18) 



i 0.070, ao = WRq, ai = 5.75 x 10" Rq = 
0.13pc and m ~ 1.2. This distribution implies that 
the numbers of wide binary systems per logarithmic in- 
terval are equal, and that approximately 50 percent of 
the stellar systems have orbital p eriods less than 100 yr 
(iHan. Podsiadlowski fc Eggleton| [l995') . 

(iv)We simply assume a single starburst (i.e. W^^ Mq 
in stars are produced one time) or a constant star formation 
rate 5* (SFR) over the last 15 Gyr calibrated so that one 
bina ry with Mi > 0.8Mq is born in the G alaxy each year 
(see Iben fc Tutukovlll984l : iHan. Podsiadlows ki & Egglcto^ 
ll995l : lHurlev et aLlbOO^ T From this calib r ation , we can get 
5* = 5 yr~^ (see also IWiUems fc Kdbl \2004 ). The con- 
st ant_staj_Jonna^ is consistent with the estimation 
of iTimmes et all l| 19971 ). which successfully reproduces the 
^®A1 1.809-MeV gamma-ray line and the core-c ollapse su- 
pernova rate in the Galaxy (|Timmes et al.|[l997h . 



5 THE RESULTS OF BINARY POPULATION 
SYNTHESIS 

5.1 the Birth rates of SNe la 

Fig. |6] shows the Galactic birth rates of SNe la (i.e. 
Z = 0.02 and SFR= 5.OAf0/yr) from the WD-hMS 
channel. From the figure, we see that the Galactic birth 
rate is around 0.7-l.Ox 10~ '^yr~^, consistent with that of 
iHan fc Podsiadlowskil l|2004l ). This result is lower but com- 
parable to that inferred from obse rv ations (3-4x 10~'^yr" 



van den Bergh fc TammannI Il99ll : ICappellaro fc Turattd 
1997f) . 



The birth rate of SNe la at various metallicities are pre- 
sented in Fig.[7l In this figure, we see that most supernovae 
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Figure 6. The evolution of the birth rates of SNe la for a constant 
star formation rate (Z=0.02, SFR=5MQyr~^). Solid and dashed 
lines show the cases with acE = 3.0 and «CE = l-Oi respectively. 




log(l/yr) 

Figure 7. The evolution of the birth rates of SNe la for a single 
starburst of 10^^ Mq for different metallicities. The upper panel 
shows the cases with ocE = 3.0, while the cases of acE = 1.0 are 
in the bottom panel. 



occur between 0.2 Gyr and 2 Gyr after the starburst. An 
interesting phenomenon is that a high metallicity leads to a 
systematically earlier explosion time of SNe la, which due 
to the effect of metallicity Z on the maximum initial mass 
of the companion for SNe la and on the stellar evolution. 
As shown in Fig. |4l M2 increases with metallicity Z. Gen- 
erally, a massive star evolves more quickly than a low-mass 
one. So, the explosion time is earlier with Z from this view. 
Though a high Z also slows the evolution of a star down, 
its influence is much less than that of st ellar mass based on 
detailed calculation s of stellar evolution l|Umeda et al ] | 19991 : 
IChen fc Tou^l2007l ). 

We also see in Fig.[7]that the peak value increases with 
metallicity Z. This comes from the fact that the parameter 



space for SNe la increases with metallicity, e.g. the range of 
the initial masses of WDs is larger for a high metallicity. In 
the upper panel of Fig. [T] however, the peak value for the 
case of ^ = 0.0001 is higher than that of Z = 0.001, which 
is mainly from the influence of acE and Z on the TPAGB 
channel (see subsection I5.2.ip . acE significantly affects the 
peak values only at low metallicities, that is, the peak value 
of the birth rate is obviously larger for acE = 1.0 in com- 
parison to QcE = 3.0 when Z < 0.004. These phenomena 
originate from the infiuences of oce and metallicity on the 
He star channel and the TPAGB channel (see subsection 
ISXTl for details). 

However, the result s abo ve are opposite to the re- 
cent study of IGuo et al. (|2008|). where an a nalytic method 
was used (s e e also iHachisu et al.l[l999al lbl and lNomoto et all 
1 19991 . 12OO3I I. IGuo et al.l lj2008 ll found that the peak of the 
birth rate for a low Z appears significantly earlier than that 
for a high Z. Also, they have not found the dependence of 
the peak value of birth rate on metallicity. This discrepancy 
comes from the different approaches in calculating IM2I in 
the two studies. As described in section [J] the value of IM2I 
is crucial for the accretion of CO WDs. Thus, initial param- 
eter spaces for SNe la are different in the two studies, i.e. 
a low metallic ity leads to a hig her maximum mass of the 
progenitors in IGuo et aP (|2008l ) . This tendency is the op- 
posite of that in this paper. Then, the delay times of SNe 
la, which are closely relevant to the initial companion mass, 
are also different. The final birth rates are therefore dif- 
ferent. For example, given Z = 0.02, the peak value from 
the analytic method is larger than that from detailed bi- 
nary evolution calc ulations by a factor of 3, if a single star- 
burst is assumed (|Han fc Podsiadlowskil |2004| ) . Note that 
our results are more physical tha n those obtaine d from the 
simple analyti c met h od us ed b y iHachisu et al.| (Il999al lbl). 
iNomoto et all (|l999l . 120031 ') and IGuo et all (|2008h . as our 
results are based on detailed binary evolution calculations 
with the latest input physics. 



5.2 Distribution of initial parameters of WD + 
MS systems for SNe la 

Observationally, some WD -I- MS syst ems are possible pro- 
genit ors of SNe la (see the review of IParthasarathv et al.l 
[2003). Further studies are necessary to finally confirm them 
(from both observations and theories). In this section, we 
will present some properties of initial WD -I- MS systems 
for SNe la , which may help us to search for the potential 
progenitors of SNe la in various environments. 



5.2.1 Distribution of initial masses of WDs 

The distribution of initial WD masses from different chan- 
nels are shown in Fig. [S] in which the relative importance 
of the three channels is clearly shown. In Fig. |9l we present 
the distribution of the initial masses of the CO WDs leading 
to SNe la for various Z and acE. Fig. [S] shows that acE is 
crucial for the He star channel. If a system experiences CE 
phase before rather than after helium burning, it will more 
likely merge due to a large binding energy and a short pri- 
mordial orbital period, especially for qce = 1.0. Thus, the 
contribution of WD + MS systems from the He star channel 
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Figure 8. The distribution of tlie initial WD masses in WD + MS systems from different SN la channels for various metallicities and 
different ciCE ■ The red, green and blue lines represent those from He star, EAGB and TPAGB channel, respectively. 



is less when acE = 1-0 and increases when qce = 3.0. Since 
the stellar radius increases with metallicity Z, a primordial 
binary with given Mf , M2 and is more likely to undergo 
He star channel evolution if it has a high Z. The importance 
of the He star channel thus increases with metallicity Z. Ac- 
tually, no WD + MS systems result from this channel when 
Z < 0.004 in our simulation. He star channel has an im- 
portant contribution to WD + MS systems with CO WDs 
around O.TSM© (see the peaks at low masses in Fig. |9]). 

The EAGB channel is not significantly infiuenced by 
acE for its moderate primordial orbital period, and it is the 
only contributor to the peaks at masses^ O.92M0 in the 
upper panel of Fig. O EAGB channel also produces many 
WD -I- MS system with CO WD mass around 0.78Mq (see 
Fig ID). 

The massive CO WDs (i.e the high-mass tail for Z > 
0.004, the plateau from 0.9Mq to I.2OM0 for Z ^ 0.004 
with acE = 1.0, and peak around I.ISM© with qce = 3.0 



and Z = 0.0001 in Fig. are mainly from TPAGB channel 
(see the blue lines in Fig. [8]). Because of the low binding 
energy of the common envelope and a long primordial or- 
bital period, acE has a remarkable influence on CO-fWD 
systems from the TPAGB channel. Generally, if a CE can 
be ejected, a low acE produces a shorter orbital-period WD 
+ MS system, which is more likely to fulfill the conditions 
for SNe la. Therefore, we see obvious contributions from 
the TPAGB channel when qce = 1-0, but no contribution 
from this channel when qce = 3.0 except for the case of 
Z = 0.0001. We explain this as follows: For a low Z, due to 
relatively small stellar radius, binaries undergoing TPAGB 
channels usually have shorter primordial orbital periods and 
tighter common envelopes, and the produced CO WD -I- MS 
binaries have shorter orbital periods after CE ejection and 
are more likely to contribute to SNe la. This effect increases 
with the decrease of metallicity. We therefore see the high- 
mass peak in the figure for metallicity as low as 0.0001 with 
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QCE = 3.0. It is a little bit different with acE = 1.0, for 
which the produced WD + MS systems have shorter initial 
orbital periods in comparison to acE = 3.0. Thus, many WD 
+ MS systems, which can not lead to SNe la with acE = 3.0, 
contribute to SNe la with acE = 1.0, and produce a plateau 
of CO WD masses instead of a peak (the bottom panel in 
Fig.®. 

As discussed above, for high Z, the contribution to 
SNe la from the He star channel increases, while that from 
the TPAGB channel decreases with acE, which leads to a 
nearly constant birth rate of SNe la for different values of 
acE. However, for the cases with low metallicities, the He 
star channel has no contribution to SNe la at either value 
adopted for acE, while a large amount of WD + MS systems 
from the TPAGB channel contribute to SNe la at the low 
acE, resulting in a higher birth rate in comparison to the 
high acE. 



I ' I ' ' I 




5.2.2 Distribution of initial secondary masses 

Fig.[lO]presents the distributions of the initial masses of sec- 
ondaries for SNe la for various metallicity Z and acE. The 
distributions for different metallicities have similar shapes 
(i.e. a low-mass sharp peak with a high-mass tail). Since the 
contours for SNe la move to higher masses with Z, the peak 
mass of the secondaries also moves to higher masses with Z. 
The difference of the peak initial masses of the secondaries 
between Z = 0.03 and Z = 0.0001 is as large as ~ 0.7 Mq. 



5.2.3 Distribution of initial orbital periods 

The distributions of the initial orbital periods are shown 
in Fig. 111! Similar to the distribution of initial masses of 
CO WDs, there are also double peaks for Z = 0.0001 and 
cecE = 3.0, which correspond to the low-mass peak and the 
high-mass one in Fig. O respectively. 



5.2.4 Distribution of initial separations 

The distributions of the initial separations are shown in Fig. 
1121 Similar to the distributions of the initial masses of CO 
WDs and that of the initial orbital periods, there are dou- 
ble peaks for the case of Z — 0.0001 and oce = 3.0. The 
low-separation peak is from EAGB channel and the high- 
separation one is from TPAGB channel. 

We see from Fig. 1121 that there are almost no initial 
systems with separations larger than ~ 26Rq . This is a nat- 
ural result since we only focus on WD -I- MS systems in 
this paper. For systems with long initial separations, mass 
transfer occurs on RGB, which is beyond the scope of this 
paper. 



6 DISCUSSION 

6.1 Correlation between the peak luminosity of 
SNe la and metallicity 

In our study, the masses of CO WDs leading to SNe la 
are more massive in a low-metallicity environment (see Figs 
[5] and [9]). Some previous studies showed that a massive 
CO WD leads to a lower C/0 ratio, and thus a lower 



Figure 9. The distribution of the initial WD masses in WD 4- MS 
systems which can ultimately produce SNe la for various metal- 
licities and different acE- 




Figure 10. The distribution of the initial masses of secondaries 
in WD + MS systems which can ultimately produce SNe la for 
different metallicities and different acE- 



amount of ^^Ni synthesize d in the thermonuclear explosion 
(|Nomoto et alJ U999. 2003 ) , which results in a lower lurninos- 
ity of SNe la (|Arnett.l98a : lArnett. Branch fc Wheeleij[l985l : 
lBranchlll992l '). According to this, our results means that the 
mean maximum luminosity of SNe la increases with metal- 
licity Z, implying that the average maximum luminosity de- 
creases with redshift, since metallicity generally decreases 
with redshift. Observationally, the mean peak brightness of 
SNe la in a galaxy has less variation in the outer regions tha n 
in the inner r egions (Wans et al.) Il997l : iRiess et al.1 Il999l ). 
iNomoto et all (|1999 , .2003 ) suggested that this results from 
the effect of metalhcity, i.e. the maximum peak brightness 
of SNe la is larger in the inner region than that in the outer 
region of a galaxy since m etallicity decreases along the ra- 
dial direction of a galaxy (|Baes et al.ll2007l : iLemasle et al.l 
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licities may be u p to 0.3 ( see F iKO, consistent with 
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Figure 11. The distribution of the initial orbital periods of the 
WD+MS systems which can ultimately produce SNe la for dif- 
ferent metallicities and different ocE- 
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Figure 12. The distribution of the initial separations of the 
WD+MS systems which can ultimately produce SNe la for dif- 
ferent metallicities and different «CB- 



[20o3), while the minimum peak brightness of SNe la is sim- 
ilar because the maximum initial CO WD mass is almost 
equal, i.e. ~ _L2Mq (see Fig. [5] in this paper and Fig. 1 in 
iMeng et al.ll200^ . which leads to the lowest C/0 ratio of 
a CO WD and thus the least amount of ^^Ni synthesized 
in the ther monuclear explosion (iNomoto et al ]Ll999, 2003). 
In addition, iHamuv et all l|l996bl ) did not find the classical 
Malmquist bias in their sample and the most distant SNe 
la are not significantly brighter than nearby samples. This 
phen omenon is at least not inconsistent with our prediction. 

iLanger et al.l ((2000,) suggested that on average, the ini- 
tial masses of CO WDs for SNe la with a low metallicity 
{Z = 0.001) is larger than that with a high metallicity 
(Z = 0.02) by about 0.2 M©. From the BPS study in this pa- 
per, the difference of the CO WD masses from various metal- 



iP „ _ , 

the suggestion of iLanger et al. ||2000|). Detailed single steL 
lar evolution calculations (*Han. Podsiadlo wski fc EggletonI 
[1994; Girar di et al. 2000: Meng ct al. 2008 ) show that, for a 
given main-sequence star, the CO WD mass decreases with 
metallicity Z, that is, a 0.2 Mq difference of CO WD masses 
can be obtained from different metallicities l|Meng et al.l 
[2008?). Since the CO WD masses from single stellar evolution 
may be taken as upper limits of CO WD masses from binary 
evolution to some extent, the mass difference of CO WDs 
for SNe la from different metallicties might be achieved from 
the evolution prior to CO WD formation. 



6.2 Other possible channels for SNe la 

From our study, the WD -I- MS channel 
can only acco unt for about 1/3 of the SNe 

la observeddvan den Bergh fc TammannI Il99ll : 

ICappellaro fc Turattol 1 19971 ). Therefore, there may be 
other channels or mechanisms contributing to SNe la. As 
mentioned in section [TJ a wide symbiotic system, WD 
+ RG , is a possible progenitor of SNe la (|Hachisu et al-l 
Il999bl ). though previous BPS studies indicated little 
contributions from this channel (|Yungelsqn fc Livi d Il998l: 
[Han & Podsiadlowski 200 4). Recently, [Hachisu fc Katd 
( 200a l2006al lbD and iHachisu et all (|2007l ) suggested that 
several recurrent novae are probable progenitors of SNe 
la and som e of them b elong to the WD-I-RG channel. 
Meanwhile, IPatat et al.l l|2007l ) suggested that the com- 
panion of the progenitor of SN 2006X may be an early 
RGB star. So, a further study of this channel is necessary. 
An alternative is the double-degenerate ( DD) channel 
(|lben fc Tutukovl Il984l : IWhelan fc lbeiill 19871). although it 
is theoretically less favored i Hillebrandt fc NiemeveJl200Gl ). 
In this channel, two CO WDs with a total mass larger than 
the Chandrasekhar mass limit may coalesce and explode as 
a SN la. The birth rate from this channel is comparable to 
the observational value ([Haij|l998l : lYungelson fc Liviolll998l , 
l200d : iTutukov fc Yungelsonll2002D . and SN 2003fg and SN 
2005hi likely resulted from the DD channel (iHowell et al.l 
l200d : lBranch|[200^ : lOuimbv. Hoflich fc Wheeled l2007l ). Ob^ 
servationally, a large am ount of DD systems are discovered 
dNapiwotzki et al.l [200j ). but only KPD 1930-1-2752 is a 
possible progenitor candidate for a SN la via DD channel 
(Geier et al., ,2007, ). The total mass of KPD 1930-1-2752 
(~ I.52M0) exceeds the Chandrasekhar mass limit and the 
time scale of coalescence is about 200 Myr estimated from 
orbital shrinkage caused b y gravitational wave radiation 
(Ccie r et al.ll2007l ). However. lErgma. Fedorova fc YungelsonI 
argued that, from detailed binary evolution calcula- 
tion, the final mass of KPD 1930-1-2752 is smaller than the 
Chandrasekhar mass limit due to a large amount of mass 
loss during evolution. In addition, earlier numerical simula- 
tions showed that the most probable fate of the coalescence 
is an accretion-induced coll apse and, finally, neutron sta r 
formation (see the review bv lHillebrandt fc Niemeverll200d ) . 
A definitive conclusion for DD model is thus premature at 
pres ent, and furth er studie s are needed. 



Liebert et al 



■ I2OO5I ) and 

IWickramasinghe fc Ferrariol ( 20051 ) found that about 
10% of WDs have magnetic fields higher than IMG. The 
mean mass of these WDs is 0.93 Mq, compared to mean 
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mass of all WDs which is 0.56 Mq (see the review by 
IParthasarathv et aLll2007l for details). Thus, the magnetic 
WDs are more likely to reach the Chandrasekhar mass limit 
by accretion. Meanwhile, the magnetic field may also affect 
some properties of WD+MS systems, e.g. the mass transfer 
rate, the critical accretion rate, the thermonuclear reaction 
rate etc, leading to a different birth rate of SNe la. 



6.3 2002ic-like supernovae 

Until SN 2002ic was discovered (|Hamuv et all 120031 '). it 
was long believed that there are no hydrogen lines in the 
spectra of SNe la. The strong hydrogen lines in the spec- 
tra of SN 2002ic were explained by the interaction be- 
tween the SN eject a and the circumstellar material (CSM) 
ijHamuv et al.ll2003l l. The CS M is aspheric around the SN la 



from spectropolarimetry data (jPeng et al.ll2004 Wang et"aLl 



20041) and has a mass of 0.5-6M( 



I Wang ^ 

Chugai fc Yungelsonll2004l : lUenishi et al.ll2004l : iKotak et all 
20041 '). The shape of the light curve showed that there 



et al 



g et al.l 
1 l2004l: 



was a delay between the explosion and the interaction, 
indicating a cavity betw e en the progenitor and the CSM 
IWood-Vasev et all |2004| : IWood-Vasev fc Sokoloskil l200d ). 
Recently, two co-twins of SN 2002ic (SN 2005 gi and SN 
2006g y ) were also found (i AIdCTing_clal: 20Q6i ; lofek et all 



I2OO7I ). iBenetti et all (HoO^argued that SN 2002ic is not 
a Type la supernova, but a Type Ic SN surrounded by a 
structured H-rich CSM, possibly asymmetric. However, re- 
cent spectral comparison between SN 2005gj, SN 2002ic and 
SNe la provided f urther evidence of Type la rather than 
Ic for SN 2002ic (|Prieto et al.l 120071 '). Thus, we still take 
SN 2002ic as a SN la. To explain these rare objects, many 
models were suggested and here we just list some of them: 
Hachisu et"ai] Il999bl: iHaniuy et"aLl l2003l: iLivio fc Riesi 



20031: IChugai fc YungelsonI |2004|: iHan fc Podsiadlowski 



200e 



Among all the models listed above, the model of 
iHan fc Podsiadlowskil (|2006l ) is the best model to match 
with the observation at present, especially about the 
birth rate of these rare objects an d the delay time 
jAldering et al.ll2006l: jPrieto et al.|[200'iD . In the scenario of 
iHan fc Podsiadlowskil 120061 ). SN 2002ic might be from the 
WD -I- MS channel, where the CO WD accretes material 
from its relatively massive companion (~ S.OMq), and in- 
creases its mass to ~ l.SOM© before experiencing a delayed 
dynamical instability. From this scenario, we examined the 
initial parameters for 2002ic-like SNe la by assuming that 
the CO WD can increase its mass to 1.378 Mq and explode 
as a SN 2002ic-like case if the mass of the CO WD exceeds 
1.30 Mq before the delayed dynamical instability. As shown 
in Figs. and ijfl. the parameter space for the SN 2002ic-like 
case becomes larger with metallicity. For Z = 0.0001, no sys- 
tems become 2002ic-like supernovae. This means that the SN 
2002ic-like case could not oc cur in extremely low-metallic ity 
environments if the model in \Han & Podsiadlowski Ij200(i ) is 
appropriate. Observationally, the host galaxy of SN 2006gy 
has a solar metallicity, while there is no information on the 
metallicity of SN 2002ic since any association between SN 



2002ic and the galaxies near SN 2002ic have been ruled out 
iHamuv et al.l l2003l : iKotak et all 120041 '). However, we may 
speculate that the metallicity of SN 2002ic is likely to be 
higher than 0.0001 from the redshift oi z = 0.066. For 
the host galaxy of S N 2005gj, we only know Z/Zq < 0.3 
(|Aldering et al.ll2006l '), which might not be as low as 0.0001. 
Our model predicts that no SN 2002ic-like events would be 
observed in extremely low metallicity environments or at 
very high redshift, which can be tested with an enlarged 
sample of 2002ic-like supernova in the future. 



7 SUMMARY AND CONCLUSION 



^ See my personal web site\lrittp://www.ynao.ac.cn/f^ bps /down- 
load/xiangcunmcng.htm to see the cases with other metallicities. 



Adopting the prescription of iHachisu et al.l (|l999al l for 
the mass accretion of CO WDs and assuming that the 
prescription is valid for all metallicities, we have studied 
the progenitors of SNe la in the single degenerate channel 
(WD+MS) via detailed binary evolution calculations. The 
initial parameters for SNe la in the (logP',M2) plane 
are obtained for ten metallicities. One can download the 
FORTRAN code for the contours leading to SNe la at 
\http://www.ynao.ac.cn/\ " bps/ download/ xiangcunmeng.htm . 

We found that the contours for SNe la move to high 
mass and long orbital period with increasing metallicity. The 
minimum initial mass of the CO WDs for SNe la sharply de- 
creases with metallicity and the difference of the minimum 
initial mass between Z = 0.06 and Z — 0.0001 is as large as 
0.24 Mq. Incorporating the binary evolution calculation re- 
sults in this paper into Hurley's rapid binary evolution code, 
we have studied the evolution of the birth rates of SNe la 
with time. The Galactic birth rate from the WD-|~MS chan- 
nel is lower than but still comparable, within a factor of a 
few, to that inferred from observations. For the cases of a 
single starburst, the SNe la explosions occur earlier and the 
peak value of the birth rate is larger for a high metallicity. 
The distributions of the initial masses of CO WDs, the ini- 
tial masses of secondaries, the initial periods and the initial 
separations e volve s ignificantly with metallicity. Based on 
the model of iHan fc Podsiadlowski (2006), we predict that 
2002ic-like supernovae would not occur in extremely low- 
metallicity environments. 



ACKNOWLEDGMENTS 

We thank Dr. Richard Pokorny for his kind help in improv- 
ing the language of this paper. This work was in part sup- 
ported by Natural Science Foundation of China under Grant 
Nos. 10433030, 10521001, 2007CB815406 and 10603013. 



REFERENCES 

Aldering G., Antilogus P., Bailey S. et al., 2006, ApJ, 650, 
510 

Alexander D. R., Ferguson J. W., 1994, ApJ, 437, 879 
Arnett W.D., 1982, ApJ, 253, 785 

Arnett, W.D., Branch, D., Wheeler, J.C., 1985, Nature, 
314, 337 

Baes M., Sil'chenko O.K., Moiseev A.V. et al., 2007, A&A, 
467, 991 



14 Meng, Chen and Han 



Benetti S., Cappellaro E., Turatto M. et al., 2006, ApJ, 

653, L129 
Branch D., 1992, ApJ, 392, 35 
Branch D., Bcrgh S.V., 1993, AJ, 105, 2231 
Branch D., 2004, Nature, 431, 1044 
Branch D., 2006, Nature, 443, 283 

Cappellaro E., Turatto M., Tsvetkov D.Y., Bartunov O.S., 

PoUas C, Evans R., Hamuy M., 1997, A&A, 322,431 
Cappellaro E., Turatto M., 1997, in Ruiz-Lapuente P., Can- 

nal R., Isern J., cds. Thermonuclear Supernovae. Kluwer, 

Dordrecht, p. 77 
Chugai N.N., Yungelson L.R., 2004, Astronomy Letters, 30, 

65 

Chen X., Han, Z., 2002, MNRAS, 335, 948 

Chen X., Han, Z., 2003, MNRAS, 341, 662 

Chen X., Tout C.A., 2007, ChJAA, 7, 2, 245 

Chen X., Han, Z., 2008, MNRAS, 387, 1416, arXiv: 

0804.2294 

Deng J., Kawabata K.S., Ohyama Y. ot al., 2004, ApJ, 605, 
L37 

Di Stcfano R., Kong A.K.H., 2003, ApJ, 592, 884 
Eggloton P.P., 1971, MNRAS, 151, 351 
Egglcton P.P., 1972, MNRAS, 156, 361 
Eggleton P.P., 1973, MNRAS, 163, 279 
Eggleton P. P., Tout C. A., Fitechett M. J., 1989, ApJ, 347, 
998 

Ergma E., Fedorova A.V., Yungelson L.R., 2001, A&A, 376, 
L9 

Guo W., Zhang F., Meng X., Li Z., Han Z., 2008, ChJAA, 
8, 63 

Gcior S., Ncsslingcr S., Hcbcr U., Przybilla N., Napiwotzki 

R., Kudritzki R.-R, 2007, A&A, 464, 299 
Girardi L., Bressan A., Bertelh G., 2000, A&A, 141, 371 
Goldberg D., Mazeh T., 1994, A&A, 282, 801 
Hachisu I., Kato M., Nomoto K., ApJ, 1996, 470, L97 
Hachisu I., Kato M., Nomoto K., Umeda H., 1999a, ApJ, 

519, 314 

Hachisu I., Kato M., Nomoto K., 1999b, ApJ, 522, 487 

Hachisu I., Kato M., 2003a, ApJ, 588, 1003 

Hachisu I., Kato M., 2003b, ApJ, 590, 445 

Hachisu I., Kato M., 2005, ApJ, 631, 1094 

Hachisu I., Kato M., 2006a, ApJ, 642, L52 

Hachisu I., Kato M., 2006b, ApJ, 651, L141 

Hachisu L, Kato M., Luna G.J.M., 2007, ApJ, 659, L153 

Hachisu I., Kato M., Nomoto K., 2008, ApJ, inpress, 

(arXiv: 0710.0319) 
Hamuy M., Phillips M.M., Schommer R.A., Schommer 

R.A.,Suntzeff N.B., Maza J., Aviles R., 1996, AJ, 112, 

2391 

Hamuy M., Phillips M.M., Suntzcff N.B, Schommer R.A. 

1996b, AJ, 112, 2398 
Hamuy M. et al., 2003, Nature, 424, 651 
Han Z., Podsiadlowski P., Eggleton P.P., 1994, MNRAS, 

270, 121 

Han Z., Podsiadlowski P., Eggleton P.P., 1995, MNRAS, 
272, 800 

Han Z., 1998, MNRAS, 296, 1019 

Han Z., Tout C.A., Eggleton P.P., 2000, MNRAS, 319, 215 
Han Z., Podsiadlowski Ph., Maxtcd P. F. L., Marsh T. R., 

Ivanova N., 2002, MNRAS, 336, 449 
Han Z., Podsiadlowski Ph., 2004, MNRAS, 350, 1301 
Han Z., Podsiadlowski Ph., 2006, MNRAS, 368, 1095 



Hillebrandt W., Niemeyer J.C., 2000, ARA&A, 38, 191 

Hjellming M.S., Webbink R.F., 1987, ApJ, 318, 794 

Howell D.A. ct al., 2006, Nature, 443, 308 

Hurley J.R., Pols O.R., Tout C.A., 2000, MNRAS, 315, 543 

Hurley J.R., Tout C.A., Pols O.R., 2002, MNRAS, 329, 897 

Iben I., Tutukov A.V., 1984, ApJS, 54, 335 

Iglesias C. A., Rogers F. J., 1996, ApJ, 464, 943 

Ihara Y., Ozaki J., Doi M. et al., 2007, PASJ, 59, 811, 

arXiv: 0706.3259 
Kato M., Hachisu L, 2004, ApJ, 613, L129 
Kippenhahn R., Weigert A., 1967, ZA, 65, 251 
Kobayashi C, Tsujimoto T., Nomoto K. et al., 1998, ApJ, 

503, L155 

Kotak R., Moiklc W.P.S., Adamson S. ot al., 2004, MN- 
RAS, 354, L13 

Langer N., Deutschmann A., Wellstein S. et al., 2000, A&A, 
362, 1046 

Lcibundgut B., 2000, A&ARv, 10, 179 

Lemasle B., Piersimoni A., PediccUi P. et al., 2007, arXiv: 

0711.3988 
Leonard D.C., 2007, ApJ, 670, 1275 
Li X.D., van den Heuvel E.P.J. , 1997, A&A, 322, L9 
Liobcrt J., Bergeron P., Holbcrg J.B., 2003, A J, 125, 348 
Liebert J., Bergeron P., Holberg J.B., 2005, ApJS, 156, 47 
Livio M., Soker N., 1988, ApJ, 329, 764 
Livio M., Riess A., 2003, ApJ, 594, L93 
Mannucci F., Delia Valle M., Panagia N., Cappellaro E., 

Crcsci G., Maiolino R., Petrosian A., Turatto M., 2005, 

A&A, 433, 807 
Mannucci F., Delia Valle M., Panagia N., 2006, MNRAS, 

370, 773 

Mattila S., Lundqvist P., SoUerman J. ct al., 2005, A&A, 
443, 649 

Mazeh T., Goldberg D., Duquennoy A., Mayor M., 1992, 
ApJ, 401, 265 

Meng X., Chen X., Tout C.A., Han Z., 2006, ChJAA, 6, 4, 
461 

Meng X., Chen X., Han Z., 2008, A&A, 487, 625, arXiv: 
0710.2397. 

Miller G.E., Scalo J.M., 1979, ApJS, 41, 513 
Napiwotzki R., Karl C, Nelemans G. et al., 2004, RMxAC, 
20, 113 

Nelemans G., Verbunt F., Yungelson L.R. et al., 2000, 

A&A, 360, 1011 
Nelemans G., Tout C.A., 2005, MNRAS, 356, 753 
Nomoto K., Thielemann F-K., Yokoi K., 1984, ApJ, 286, 

644 

Nomoto K., Umeda H., Hachisu L Kato M., Kobayashi C, 

Tsujimoto T., 1999, in Truran J., Nicmoycr T., cds, Type 

la Suppcrnova :Theory and Cosmology. Cambridge Univ. 

Press, New York, p. 63 
Nomoto K., Uenishi T., Kobayashi C. Umeda H., Ohkubo 

T., Hachisu L, Kato M., 2003, in Hillebrandt W., Lei- 

bundgut B., eds, From Twilight to Highlight: The Physics 

of supernova, ESO/Spririgcr serious "ESQ Astrophysics 

Symposia" Berlin: Springer, p. 115 
Ofek E.G., Cameron P.B., Kaslwal M.M. et al., 2007, ApJ, 

659, L13, arXiv: 0612408 
Paczyriski B., 1976, in Egglcton P.P., Mitton S., Whclan J., 

eds. Structure and Evolution of Close Binaries. Kluwer, 

Dordrecht, p. 75 
Parthasarathy M., Branch D., Jeffery D.J., Baron E., 2007, 



The progenitor of type la supernovae with different metallicities 15 



NewAR, 51, 524, arXiv: 0703415 
Patat E. et al.. Science, 317, 924 

Podsiadlowski P., Rappaport S., Pfahl, 2002, ApJ, 565, 
1107 

Podsiadlowski P., Mazzali P. A., Lesaffre P., Wolf C, 

Forster F., 2006, arXiv: 0608324 
Perlmutter S. et al., 1999, ApJ, 517, 565 
Phillips M.M., 1993, ApJ, 413, L105 

Pols O.R., Tout C.A., Eggleton P.P. et al., 1995, MNRAS, 
274, 964 

Pols O.R., Tout C.A., Schroder K.P. et al., 1997, MNRAS, 
289, 869 

Pols O.R., Schroder K.P., Hurly J.R. et al., 1998, MNRAS, 

298, 525 

Pricto J.L. et al., 2008, ApJ, 673, 999, arXiv: 0707.0690 
Prieto J.L. et al., 2007, arXiv: 0706.4088 
Quimby R., P. Hoflich, J.C. Wheeler, 2007, ApJ, 666, 1083 
Riess A. et al., 1999, AJ, 117, 707 
Riess A. et al., 1998, AJ, 116, 1009 
Ruiz-Lapuente P. et al., 2004, Nature, 431, 1069 
Scannapieco E., Bildsten L., 2005, ApJ, 629, L85 
Schroder K.P., Pols O.R., Eggleton P.P., 1997, MNRAS, 
285, 696 

Shanks T., Allen P.D., Hoyle F. et al., 2002, ASPC, 283, 
274 

Timmes F.X., Diehl R., Hartmann D.H., 1997, ApJ, 479, 
760 

Timmes F.X., Brown E.F., Truran J.W., 2003, ApJ, 590, 
L83 

Travaglio C, Hillebrandt W., Reinecke M., 2005, A&A, 
443, 1007 

Tutukov A.V., Yungclson L.R., 2002, Astron. Rep., 46, 667 
Totani T., Morokuma T., Oda T. et al., 2008, PASJ, 60, 

1327, arXiv: 0804.0909 
Uenishi T., Suzuki T., Nomoto K. et al., 2004, Rev. Mex. 

A&A, 20, 219 

Umcda H., Nomoto K., Yamaoka H. ct al., 1999, ApJ, 513, 
861 

van den Bergh S., Tammann G.A., 1991, ARA&A, 29, 363 
Voss R., & Nelemans G., 2008, Nature, 451, 802 
Wang L., Hoflich P., Wheeler J.C, 1997, ApJ, 483, L29 
Wang L., Baadc D., Hoflich P. et al., 2004, ApJ, 604, L53 
Webbink R. F., 1988, in The Symbiotic Phenomenon, cds. 

J. Mikolajewska, M. Friedjung, S. J. Kenyon & R. Viotti 

(Kluwer: Dordrecht), p. 311 
Webbink R., 2007, arXiv: 0704.0280 
Whelan J., Iben 1., 1973, ApJ, 186, 1007 
Whelan J., Iben I., 1987, in Philipp A.G.D., Hayes D.S., 

Liebert J.W., eds, lAU Colloq.95, Second Conference on 

Faint Blue Stars. Davis Press, Schenectady, p. 445 
Wickramasinghe D.T., Ferrario L., 2005, MNRAS, 356, 

1576 

Willems B., Kolb U., 2004, A&A, 419, 1057 
Wood-Vasey W.M., Wang L., Aldering G., 2004, ApJ, 616, 
339 

Wood-Vascy W.M., Sokoloski J.L., 2006, ApJ, 645, L53 
Yungelson L., Livio M., Tutukou A. Kenyon S.J., 1995, 

ApJ, 447, 656 
Yungclson L., Livio M., 1998, ApJ, 497, 168 
Yungelson L., Livio M., 2000, ApJ, 528, 108 



